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Employing the method of mapping the spin problem onto a particle one, we have derived the particle 
Hamiltonian for a biaxial spin system with a transverse or longitudinal magnetic field. Using the 
Hamiltonian and introducing the parameter p{= (f/max — i?)/(f/max — Umin)) where f/max (Umin) 
corresponds to the top (bottom) of the potential and E is the energy of the particle, we have studied 
the first- or second-order transition around the crossover temperature between thermal and quantum 
regimes for the escape rate, depending on the anisotropy constant and the external magnetic field. It 
is shown that the phase boundary separating the first- and second-order transition and its crossover 
temperature are greatly infiuenced by the transverse anisotropy constant as well as the transverse or 
longitudinal magnetic field. 
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In recent years, quantum-classical escape rate transi- 
tion in the spin system has emerged as good cajadidates 
to display first- or second-order transition(FST).El □ Such 
a system is a single domain ferromagnetic particle with 
the magnetization M whose direction is subject to the 
magnetocrystalline anisotropy. In this situation the di- 
rection M has at least two equivalent stable orientations 
separated by an energy barrier U. Even though M is 
initially directed along one of these equivalent orienta- 
tions, M can be changed by thermal activation whose 
rate is proportional to exp(— f//fcsT) at high tempera- 
ture, and by quantum tunneling at a temperature low 
enough to neglect the thermal activation. In general, 
since the tunneling rate is dominated by exp(—U/huj) 
where u characterizes the width of the parabolic top of 
the barrier hindering the tunneling process, there is a 
crossover temperature Tq from the thermally activated 
to the quantum tunneling process. Whether the transi- 
ton about the crossover temperature is first- or second- 
order one is determined by the external magnetic field 
and the magnetic anisotropy constant. By controlling 
the field and choosing the anisotropy constant in an ap- 
propriate way, we expect that there exists the crossover 
temperature Tq at the phase boundary between first- 
and second-order transition. 

Theoretical studies of-jthe transition have been around, 
for some time. AfSeckcl and Larkin and OvchinnikovO 
demonstrated that a second order transition from ther- 
mal to quantum regimes can occur at the crossover 
temperature by usi«g the standard instanton technique. 
Later, ChudnovskyQ discussed the criterion to determine 
FST based on the behavior of the period of oscillations in 
the inverted potential. Since then, based ppon the map- 
ping of a spin problem onto a particle oneja Chudnovsky, 



Garanin and MartinesETa suggested the spin system with 
the uniaxial crystal symmetry which shows FST in the 
presence of a transverse and longitudinal field. A biaxial 
spin model without jan applied field has been considered 
by J.-Q. Liang et als who demonstrated that FST is de- 
termined from the ratio of the transverse to the longitu- 
dinal anisotropy constant by using the periodic instanton 
method. Even though they presented the analytical re- 
sults without field, their approach could not be simply 
extended to the situation in the presence of field.EI Very 
recently, in an effort to treat FST of the biaxial spin syjai 
tem with a longitudinal field, Garanin and Chudnovskyll3 
used a perturbation approach and obtained the phase 
boundary between the first- and the second-order transi- 
tion which is numerically corrected by the Liang et aZ.'s 
exact value in the absence of the field. Thus, relevant 
approach to treat a biaxial spin system with a transverse 
or longitudinal field has been highly required for FST. In 
fact, complete analytical solution of the problem seems 
to be considerablely important for the possibility of FST 
in molecular magnetic systems as well as a single domain 
ferromagnetic particle with many spins. In this paper, 
employing the method of mapping of a spin problem onto 
a particle one, we obtain the particle Hamiltonian of the 
biaxial spin problem with a transverse or longitudinal 
field, wfaiyh has not been obtainable from the previous 
studics.HCl Actually, such a mapping is not a regular pro- 
cedure and its form strongly depends on the form of the 
spin Hamiltonian. Using the Hamiltonian and introduc- 
ing the dimensionless energy variable p, we will study 
the first- or second-order transition around the crossover 
temperature between thermal and quantum regimes for 
the escape rate, depending on the anisotropy constant 
and the external magnetic field, and present the ana- 



lytic form of the phase boundary separating first- from 
second-order transition and the crossover temperature at 
the phase boundary. 

Consider the biaxial spin model in a transverse field 
Hx described by the Hamiltonian 



n = ~K\\Sl 



K±S — HxSx, 



(1) 



where i^|| and K± are the longitudinal and transverse 
anisotropy constants, respectively. JChis Hamiltonian can 
be mapped onto a particle problemH which describes the 
exact correspondence between the spin wave function 



V' 



s 

M=-S 



[s + My.{s - Aiyy^^ 



i2sy. 



\SM), (2) 



where \SM) are the eigenstates of Sz, and the particle 



wave function ^!{x) 
where 



exp(-7(a;)) J2m=-s "-m exp(Af x) 
d-f{x) _ Sh^l - fc)sinh(a;) + (5 - l)fcsinh(2a;) 



dx 



1 + A:cosh(2a;) 



(3) 



5 = 5+ 1/2, k = kt/{2 + kt), kt = K^/K\\ and K 
Hx/{2SK\\). The particle Hamiltonian is 



H = - 



where l/m{x) 
S'^K\\u{x), and 

u{x) 



2m{x) dx"^ 
K\\{2 + kt)[l 



1 



U{x), (4) 

-f fccosh(2x)], U{x) = 

-{hliiinh^{x){l-kf 



(l-fc)[l + fccosh(2x)) 
—2hx{\ — k) cosh(a;) — k[2hx{\ — k) cosh(x) 
-1 + cosh(2a;) - (5cosh(2a;) - l)/(452)] 
-t-fc2[cosh(2.x) - 1 



-H(cosh2(2a;) ~ cosh(2a;) + 4)/(452)]}. 



(5) 



Here we notice that the potential and the mass without 
transverse anisotropy are reduced to the results in Rcf.El. 
Even though 5 appears in the potential, the terms as- 
sociated with 5 can be neglected in the large 5 limit 
including 5 = 10. So, we will not includej-the terms with 
1/(45'^) in the subsequent consideration. E£l 

In the quasiclassical approximation the transition rate 
becomes 



r ^ J dEW{E)eM-iE-Ur^in)/T], 



(6) 



where W{E) is the probability of tunneling at an energy 
E. Since this is defined via the imaginary-time action 
W{E) ~ e~'^'^-^\ the transition rate is approximately 
given by 



r~exp(-F^i„/r), 



(7) 



where Fmin is the minimum of the effective "free energy" 
F = E + TS{E) - C/min with respect to E. Then, writing 



F{E)/T = J^{)l^dT[{m{x)/2){dx/dT)^ + U{x) 



1/^ 
the imaginary time action is given b 

fXi{E) 






Er, 



S{E) = 2 / dx>/2Mx)\/U(x)^^, 



(8) 



where ±xi {E) are the turning points for the particle with 
energy —E in the inverted potential —U{x). We note 
that the mass is coordinate dependent, which is crucial 
in changing the boundary between first- and second-order 
transition, as will be seen. In order to determine FST in 
the crossover regime, we need to consider the behavior 
of S{E) near the top of the barrier. Since the potential 
in Eq. (ph is even, we expand the integrand in Eq. m) 
near a; = which corresponds to the ton-|p|f the barrier. 
Introducing dimensionless energy variableEla p(= (/7i„ax— 

£^)/(C^max-t/min)) whcrC C/max (Umin) COrrCSpOuds tO the 

top (bottom) of the potential, the action becomes near 
the top of the barrier 



SiE) 



ttAU 



^/TMk^^+kT) 



[p + /3p2+0(/)], 



(9) 



where U2 = -S^K\\u'^'^^{Q)/2, C/4 = S^K\\u'^^\Q)/4:\, and 



^-im-) 



1--, 

3 Vl + fc 






(10) 



Here AC/ = C/max — C^min and the second term in the 
bracket comes from the coordinate dependence of the 
mass. Using the analogy with the Landau theory of 
phase transitions and the general conditions for first- and 
second-order quantunwilassipal transition of the escape 
rate discussed in Refs.El andQ, the factor in front of p^ in 
Eq. (y) determines the boundary between the first- and 
second-order transition. 

For the model without transverse field we have U2 = 
2S'^K\\k/{l + k) and C/4 = -^2S'^K\\k{l - 5fc)/[3(l -f kf] 
from Eq. (||) . If the mass does not depend on the coor- 
dinate, the sign of the factor of C/4 determines whether 
the system becomes the first- or the second-order transi- 
tion due to the anharmonicity of the potential. However, 
as has been already noticed from Eq. (Q), the mass is 
a function of x. So, we cannot simply obtain the phase 
boundary from the anharmonicity of the potential near 

S'^Ku, the 



the top of the barrier. Now, using AC/ 
action (g) is given by 



Ml 



S{E) = {^S/^/h)[p + Pp' + 0{p')], 



(11) 



where (3 — {1~ l/kt)/8. Thus the critical value of kt is 1 
implying that smaller values of kt lead to the first-ocder 
transition which is consistent with the result in Ref.QO 
In case of the second-order transition the crossover oc- 
curs at temperature Tq = LjQ/{2n) where Wq is oscil- 
lation frequency near the bottom of the inverted poten- 
tial. To estimate the frequency, we set up the Euclidean 
Euler-Lagrange equation m(a;)x-|-m'(a;)a;^/2 — C/'(a;) = 0, 
and insert x = Xb + Sx into the equation, where xt is x- 
coordinate of the barrier. Expanding to second order in 



Sx, we have Sx + ujqSx = 0, where ujq — ■\/|C/"(0)|/m(0). 

(2) 

This yields the crossover temperature given by Tq = 
S^/K\\Kj_/tt, which leads to the crossover temperature 

at the phase boundary, T(j''^ = SK\\{2 + fct)/(37r). 

For the first-order transition the approximate form of 
the crossover temperature for small k can be analytically 
calculated from the relation Tq^' ~ AU/S{E^in)- At the 
bottom of the barrier, S'(£',„inj can be calculated directly 
from the integral expression (0) or following fromQ 



9k — f 
S{E) = iS^^—^[JC{q) - (1 - a'Ma\q)], (12) 

which is derived from the Eqs. (|g) and (|g). Here £ = 
E/S'^K, and K. and 11 arc complete elliptic integrals of 
the first and third kind with q^ = (zi — l)/(zi + 1), a^ = 
g2(l-A;)/(l + fc),andzi = [k{k-l)+£]/[k{k-l)- k£]. 

This yields T^^'^ ~ {SK\\/2)/ In [(1 + ^/^~W)/k\ , which 
is approximate form for the first-order transition in the 
region of small k. 

We now consider the model with a transverse magnetic 
field. Expanding the potential (^ in powers of x, we 
obtain 



f/2 - [S^K\\/{1 + kt)]{K + h){l - K), 



(13) 



Ui = [^2/^11/(12(1 + kf)][A{l - 5fc)(l - k)hl 

+ (-1 + 38A: - b7k^)h^ - Sk{l - bk)], (14) 

and AC/ = S'^K\\{1 — hxf ■ It is noted that the condition 
in which the barrier does not vanish is U2 > 0, i.e., hx < 
1. Then, the action near the top of the barrier becomes 
from Eq. (ffl 



SiE) = 7r5 



5(1-M'/' 



ih + h^y/^ 



[p + pp^ + Oip% 



where 



/? = 



(1 



S{h, + kty 



(hx - h+)ihx - h^), 



(15) 



(16) 



^^±(^.j _ 1 - 14fc, + k^ ± (1 + h)y/l + 34kt + kf ^^^^ 



As is shown in Fig. H, the behavior of h^ shows that /3 
is negative for fcj < 1 and h^ < h^ which corresponds 
to the first-order transition. As is shown in Fig. ||, in 
the absence of the transverse anisotropy, h^ = 1/4 is 
the critical value for FST,EJ and in the presence of very 
small transverse anisotropy we have /i+ ~ (1 -I- 7A;t/2)/4 
in which the boundary becomes wider for < kt < 0.2. 
As kt continues to increase, the region where a first-order 
transition occurs, is smaller for the biaxial model than for 
the uniaxial model. This is intuitively understood that, 
since the transverse anisotropy drives the decay of the 
metastable state, it plays the role of the transverse field in 
the uniaxial case and so, for a given small transverse field 
the region for the first-order transition decreases as the 



transverse anisotropy increases. It is also noted that /i+ 
decreases linearly for kt ^ 1, i.e., just as /i+ ~ (1 — fct)/3. 
Continuing in the present case as without transverse 
case, we have the crossover temperature for the second- 
order transition 



T,Y\kt,hx) = iSK\\/7r)^{kt + /i,)(l - K). (18) 

Using Eq. ( [l7| ) for ft,+ , the crossover temperature at the 
phase boundary between the first- and second-order tran- 
sition is written as 



f^ ^ {V3/{2tt))SK\\{1 + kt)^ht{kt)/{l-kt), (19) 



We note that 



T, 



which is illustrated in Fig. 

T^"^ ~ 5/i:||[V3/(47r)](l -f 3fct) for small h and T^ 

[SK\\/{2T:)]{l + kt) for A:< 1. 

In case of the first-order transition the crossover tem- 
perature as a function of kt for small hx is approxi- 
mately estimated from the ground-state tunneling cxpp=. 
nent given by Eq. (ph or the direct integral cxpressionEj 
and A[/ considered above, which gives 



T^^\kt.hx) 



SK\\{l~hxfl[2g{kt,hx% 



(20) 



where 



9{h,hx) 



In 



x/m^+v/T~^ 



2h, 



X arctan ( ^fct/(l + kt)\/l/h1 — 



(21) 



Simple analysis for the crossover temperature shows that 
there exists a maximum of the crossover temperature 
jjmax jj^ ^j^g regime of the second-order transition, which 



from Eq. (|18|) gives T^'^'' = [SK\\/ {2t:)]{1 + h) at 
hx = {I — kt)/2. It is noted that this is the asymtotic 
form of Tq' for /ct < 1 discussed previously. As the 
transverse anisotropy increases, /i™^'^ for the maximal 
crossover temperature decreases, (Fig. g) while T^'^^ in- 
creases linearly. (Fig. y) As is summarized in Table | and 
shown in Fig. pi as kt increases, the difference between 
jjmax g^j^j jjc decreases and becomes zero at fcj = 1 which 
is the critical value in the fieldless case. 

In the presence of a longitudinal field of the spin model 
H = — ii'||S'^ + K±Sy — HzSz, we can proceed the con- 
sideration, similarly and so, we will briefiy discuss the 
essential points. In order to obtain the relation between 
hz[= Hz/{2SK^\)] and kt at the phase boundary, we need 
to have the coefficient of p^ in the action, i.e., (3 includ- 
ing C/3 and m' (xb) as well as U2 and t/4 where Xb is the 
position of the barrier, because the potential given by 



u{x) 



1 



-Ahli^-kf 



{l-k)[l + kcosh.{2x)^ 

- k[2hz{l - k) sinh(2x) - 1 + cosh(2a;) 

- (5cosh(2x)-l)/(452)] 
+ A;2[cosh(2x) - 1 

+ (cosh2(2a;) - cosh(2x) + 4)/(4S'2)]}, (22) 



is not an even function. Following the procedure dis- 
cussed previously, we have the boundary between first- 
and second-order transition 



h = {l-hl)/{l + 2hl). 



(23) 



The ratio of two anisotropy constans, kt decreases 
parabolically for /i^ ^ 1, as fcj ~ 1 — S/i^, and linearly for 
hz ^ 1, as kt — (2/3)(l — hz). This can be understood 
from the fact that, since the height of barrier decreases as 
hz increases, the first-order transition is expected for the 
larger width of the barrier which comes from the smaller 
value of the trasverse anisotropy. The crossover tempera- 
ture Tq at the phase boundary can be obtained by using 

Eq. (11) and T^^''\kuhz) = {l/2n)^\U"{xb)\/m{xi,) in 
the second order transition, which leads to 



T, 



(c) 



{SK\\/^){i-hl)/.JTT2hl. 



(24) 



Simple analysis shows that our results are consistent with 
the ones deduced from the correptjon of the perturba- 
tive calculation performed in Rcf.E3 up to the numerical 
factors. Even though the perturbative approach is less 
justified at large value of the transverse anisotropy con- 
stant, its boundary and crossover temperature are strik- 
ingly the same trend as our analytical results obtained 
from the quasiclassical method. 

In this paper we have considered the quantum-classical 
escape rate transition of a biaxial spin system in the pres- 
ence of a transverse or longitudinal field by using the 
particle Hamiltonian mapped from the spin system. The 
coordinate dependence of the particle mass was crucial 
in changing the boundary between the first- and second- 
order transition and its boundary was greatly influenced 
by the transverse anisotropy constant and external field, 
whose effect is expected to obsecsie in future experiment 
including Fcg molecular magnet Jij. 

I am indebted to E. M. Chudnovsky for many useful 
discussions. This work was supported by the Basic Sci- 
ence Research Institute Program, Ministry of Education, 
Project No. BSRI-98-2415. 
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FIG. 2. Boundary between the first- and the second-order 
transitions, where (a) h^{kt), (b) h^'^^{kt) which gives the 
maximal crossover temperature in the quantum-classical tran- 
sition, and (c) the boundary whether the barrier vanishes or 
not. 



FIG. 3. Dependence of (a) the crossover temperature Tq 
at the phase boundary and (b) the maximum crossover tem- 
perature T^^^ on the scaled anisotropy constant kt where 
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